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[l] Targeted experiments with a comprehensive chemistry- 
climate model are used to demonstrate that seasonality and 
the location of the peak wanning of sea surface temperatures 
dictate the response of stratospheric water vapor to El Nino. 
In boreal spring, El Nino events in which sea surface temper- 
ature anomalies peak in the eastern Pacific lead to a warming 
at the tropopause above the warm pool region, and subse- 
quently to more stratospheric water vapor (consistent with 
previous work). However, in fall and in early winter, and also 
during El Nino events in which the sea surface temperature 
anomaly is found mainly in the central Pacific, the response 
is qualitatively different: temperature changes in the warm 
pool region and specifically over the cold point region are 
nonuniform, and less water vapor enters the stratosphere. 
The difference in water vapor in the lower stratosphere 
between the two variants of El Nino approaches 0.3 ppmv, 
while the difference between the winter and spring responses 
exceeds 0.5 ppmv. Citation: Garfinkel, C. I., M. M. Hurwitz, L. 
D. Oman, and D. W. Waugh (2013), Contrasting effects of Central 
Pacific and Eastern Pacific El Nino on stratospheric water vapor, 
Geophys. Res. Lett., 40, 4115-4120, doi:10.1002/grl.50677. 

1. Introduction 

[ 2 ] A comprehensive understanding of interannual vari- 
ability in the amount of water vapor entering the strato- 
sphere remains elusive. In the time -mean boreal wintertime, 
enhanced convection and diabatic heating over the Indo- 
Pacific warm pool (hereafter warm pool) region leads 
to lower temperatures near the tropopause in this region 
[Highwood and Hoskins, 1998; Fueglistaler et al., 2009]. 
These low temperatures govern the amount of water vapor 
that can reach higher in the stratosphere [Mote et al., 
1996; Fueglistaler et al., 2004; Fueglistaler and Haynes, 
2005; Oman et al., 2008]. Temperatures near the cold- 
est region at the tropopause (hereafter CPR, for cold 
point region) change on interannual timescales, however. 
As the amount of water vapor that enters the strato- 
sphere is important for stratospheric chemistry [Solomon 
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et al, 1986] and radiative balance [Forster and Shine, 
1999; Solomon et al., 2010], it is important to understand 
the factors that control its entry into the stratosphere on 
all timescales. 

[ 3 ] The El Nino-Southern Oscillation (ENSO) is the 
dominant mode of interannual variability in the Tropics, and 
it has been connected to variability of the CPR [Kiladis et al., 
2001; Scaife et al., 2003; Calvo Fernandez et al., 2004; 
Garcia et al, 2007; Xie et al., 2012; Scherllin-Pirscher et 
al., 2012, among others]. While these studies differ in the 
data used and the time period chosen, they conclude that El 
Nino (EN) leads to a Rossby wave response whereby anoma- 
lously warm temperatures are present over the warm pool 
region. This subsequently leads to more water vapor entering 
the stratosphere [Zhou et al. , 2004; Fueglistaler and Haynes, 
2005; Liang et al., 2011]. 

[ 4 ] However, no two EN events are identical, and these 
different variants (or “flavors”) of EN may lead to differ- 
ent atmospheric teleconnections in the tropical tropopause 
layer, hereafter TTL. During “traditional” EN events, warm 
sea surface temperature anomalies (SSTa) are present in the 
equatorial East Pacific. This variant of EN will be referred 
to as East Pacific warming, or EPW, in the rest of this 
manuscript. During many recent El Nino events, however, 
warm SSTa were concentrated in the central Pacific (CPW, 
or central Pacific warming) [Trenberth and Stepaniak, 2001; 
Larkin and Harrison, 2005; Yu and Kao, 2007; Ashok et al., 
2007; Kug et al., 2009; Kao and Yu, 2009]. It is important 
to understand how these two forms of EN affect CPR tem- 
peratures and stratospheric water vapor. Xie et al. [2012] 
considered this issue but focused on the tropopause level 
and below. 

[ 5 ] This paper will use modeling experiments to demon- 
strate that the moistening of the stratosphere due to EN 
occurs mainly in the late boreal spring while the EN event 
is decaying, even though the strongest impact on tropopause 
temperatures is several months earlier. In addition, they 
demonstrate that while the two variants of EN have generally 
similar imprints on the tropical tropopause temperatures near 
the warm pool, the subtle differences between them lead to 
qualitatively different impacts on stratospheric water vapor. 

2. Methods and Models 

[6] Modeling experiments with the Goddard Earth 
Observing System Chemistry-Climate Model, Version 2 
(GEOSCCM) ( Rienecker et al. [2008] and section 2.2 of 
Hurwitz et al. [2010]) are used to isolate the impact of 
CPW and EPW on tropopause temperature and water vapor. 
The GEOSCCM couples the GEOS-5 atmospheric global 
climate model (GCM) [Rienecker et al., 2008] with a com- 
prehensive stratospheric chemistry module [Pawson et al., 
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2008]. The model has 72 vertical layers, with a model top at 
0.01 hPa, and all simulations discussed here were performed 
at 2° latitude x 2.5° longitude horizontal resolution. The 
convection scheme used in GEOSCCM is based on Relaxed 
Arakawa-Schubert [Moorthi and Suarez, 1992; Rienecker 
et al., 2008; Molod et al., 2012]. The model vertical lev- 
els between 150 hPa and 50 hPa are located at 139.1 hPa, 
118.3 hPa, 100.5 hPa, 85.4 hPa, 72.6 hPa, 61.5 hPa, and 
52.0 hPa; output is plotted at standard pressure levels. The 
tropopause height is defined using a blended approach which 
takes the lower of a potential vorticity-based tropopause and 
a thermally based tropopause. 

[ 7 ] On plots of 100 hPa temperature, the CPR during EN 
is enclosed with a green contour; it bounds the coldest 10% 
of the Tropics between 30°S and 30°N. The climatologi- 
cal CPR is enclosed by a cyan line; it encloses the coldest 
10% of the tropics from 30°S to 30°N for each season and 
corresponds to 189.70 K for November-December-January- 
February (NDJF) and 190.78 K for April-May-June (AMJ). 
The 10% threshold is chosen to match the approximate 
size of the CPR as shown by Fueglistaler et al. [2004] and 
in Figures 1 and 4 of Oman et al. [2008]. Flowever, this 
is just an approximate measure of coldest temperatures at 
tropical tropopause, and results are not sensitive to choos- 
ing, e.g., the 6% or 14% threshold. Temperature anomalies 
at 85 hPa resemble quantitatively those at 100 hPa, and 
we therefore show 100 hPa anomalies only for brevity. 
Stratospheric Processes and their Role in Climate- 
Chemistry-Climate Validation Activity ( SPARC-CCMVal ) 
[2010] and Gettelman et al. [2010] grade highly the rep- 
resentation of stratospheric water vapor and lower strato- 
spheric temperatures by an earlier version of GEOSCCM as 
compared to observations and to the multimodel mean of an 
ensemble of CCMs. 

[8] In order to isolate the impact of changing SSTs, other 
external boundary conditions are held fixed. Greenhouse 
gas and ozone-depleting substance concentrations are spec- 
ified in the boundary layer to represent the year 2005 in all 
experiments discussed here. (Within the atmosphere, their 
concentrations depend on transport and chemistry and subse- 
quently feed back to the radiative calculations.) In addition, 
variability related to the solar cycle and volcanic aerosols are 
not considered. 

[ 9 ] The first group of model integrations are those used in 
Hurwitz et al. [2011] and Garfinkel et al. [2012] to isolate 
the impact of CPW and EPW on the Pacific-North America 
region and the boreal winter stratosphere. Briefly, Garfinkel 
et al. [2012] examined four 50 year time-slice simulations 
forced by repeating annual cycles of sea surface tempera- 
tures and sea ice that represent CPW, EPW, and La Nina (in 
which negative SSTa are present in the eastern and central 
equatorial Pacific; hereafter LN) events. The SSTa used to 
force the simulations are shown in Figure 2 of Garfinkel et 
al. [2012] and are repeated here in Figure 1. The CPW SSTa 
peaks in the Central Pacific while the EPW SSTa peaks in the 
Eastern Pacific. We have also examined the response in the 
CPW ideal experiment from Garfinkel etal. [2012] (in which 
the SST anomalies are identically zero outside of the tropical 
Central Pacific), and the anomalies are quantitatively similar 
to those in the CPW experiment (e.g., the pattern correla- 
tion with the water vapor anomalies shown in Figure Id is 
0.97). The key point is that these model integrations pro- 
vide many samples of the atmospheric response to identical 


SSTa and are long enough to achieve statistical robustness 
(as discussed in Garfinkel et al. [2012, section 5]). Statistical 
significance is detennined by a two-tailed Student’s t test. 

[ 10 ] The second group of model integrations is a six- 
member ensemble of GEOSCCM in which all sources of 
externally forced variability are removed except for one: the 
sea surface temperatures. The SSTs follow those observed 
in the years 1980 to 2006 [Rayner et al., 2003]. Each of 
the six differ slightly in their physical parameterizations, 
whereby two of the ensemble members include an inter- 
nally generated quasi-biennial oscillation and four do not. 
All results shown here are robust to excluding those ensem- 
ble members with a quasi-biennial oscillation. The key point 
is that these simulations isolate the impact of SST varia- 
tions from 1980 to 2006 on the TTL. A composite of eastern 
Pacific warming events is formed of winters in which the 
anomalous SSTs in the Nino3 region (150°W to 90°W, 
5°S to 5°N) exceeds its climatological value by 0.85 K. 
Results are similar if a threshold of 0.7 K or 1 K is used. 
The EPW winters included are: 82/83, 87/88, 91/92, 97/98, 
and 02/03. A similar composite of Central Pacific warming 
events is formed of winters in which the index 1.5*Nino 
4-0.5*Nino3 exceeds its climatological value by 0.5 K and 
which are not already characterized as EPW events (the 
Nino4 region is 160°E to 150°W, 5°S to 5°N). The CPW 
winters included are: 90/91, 94/95, and 04/05. Results are 
qualitatively similar if a threshold of 0.35 K is used, or if we 
use one of the other CPW definitions considered in Garfinkel 
et al. [2012] (unlike the polar stratospheric response, which 
is highly sensitive to CPW definition). Results are also not 
sensitive to forming composites based on tropical precipita- 
tion anomalies, as opposed to SST anomalies (see supporting 
information). The La Nina years used are 84/85, 88/89, 
98/99, 99/00, and 05/06; we find little sensitivity in the LN 
years composited when we use either the CP or the EP defi- 
nition. The impact of SST trends on near-tropopause trends 
in these integrations is considered in (C. I. Garfinkel et al., 
Upper tropospheric and lower stratospheric zonally asym- 
metric tropical temperature trends forced by sea surface 
temperature trends: Implications for water vapor and ozone, 
submitted to J. Geophys. Res, 2013). Nonlinearity between 
the EN and LN responses is discussed in the supporting 
information. Statistical robustness is determined by evaluat- 
ing the fraction of the ensemble members who agree on the 
sign of the anomaly. Light red or blue is added when five 
out of six of the ensemble members agree on the sign of the 
change, and dark red or blue is added when all six ensemble 
members agree on the sign of the anomaly. 

[ 11 ] Combined, the two groups of models allow a more 
comprehensive assessment of the impact of ENSO on the 
tropopause region and the lower stratosphere. The six- 
member ensemble simulates the complete life cycle of 
specific, observed ENSO events. The 50 year perpetual 
ENSO integrations can isolate the response to a given SST 
distribution with more confidence than in the six-member 
ensemble. 

3. Results 

[ 12 ] We begin with the tropical water vapor response in 
the 50 year perpetual ENSO integrations. The seasonal evo- 
lution of water vapor for the EPW experiment is shown 
in Figure lc, and the seasonal evolution for the CPW 
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Figure 1 . Sea surface temperature, tropical water vapor (20°S-20°N), and cold point temperature anomalies in the perpet- 
ual CPW and EPW GEOSCCM integrations, (a, b) Contours are shown at ±0.4, ±0.8, ±1.2, ±2, and ±3; region with large 
warming (cooling) are colored red (blue), (c-d) Contours are shown at ±0.1, ±0.2, ±0.4, ±0.7, ±1 ppmv, and regions 
with anomalies significant at the 90% (99%) level are colored orange (red) or light blue (dark blue). A solid black contour 
indicates the climatological position of the tropopause. The pattern correlation between the CPW and EPW anomalies are 
shown in the title of Figure la-ld. (e) The temperature anomaly associated with EN averaged in the cold point regions at 
100 hPa (as defined in section 2) for each calendar month. 


experiment is shown in Figure Id. Tropospheric water vapor 
is increased during the mature phase (November-February) 
of all variants of EN. However, the stratospheric water 
vapor concentrations depend crucially on the variant of EN. 
Namely, the stratosphere is drier during CPW than during 
EPW. In the annual average, the water vapor anomaly at 
70 hPa in the EPW experiment is -0.13 ppmv; in the CPW 
experiment, it is -0.33 ppmv. The difference in response 
between the CPW experiment and both the LN experiment 
and the EPW experiment is statistically significant at the 
99% level. The difference between EPW and CPW peaks 
at 0.3 ppmv for certain calendar months. Furthermore, pro- 
nounced seasonality is evident in the water vapor response. 
For both the EPW and CPW experiments, moister air parcels 
enter the stratosphere during spring (April- June) than during 
other calendar months. Even though ENSO events typically 
peak in early winter, the moistening of the lower stratosphere 
noted by previous work [e.g., Zhou et al., 2004; Fueglistaler 
and Haynes, 2005; Liang et al., 2011] is delayed until late 
spring. In the rest of this paper, we will address two ques- 
tion: (1) Why is the plume of increased water vapor that is 
apparently caused by both EN variants delayed until spring? 
and (2) Why is there more stratospheric water vapor during 
EPW than during CPW? 

[n] In the interest of brevity, we will answer these ques- 
tions by focusing on the response in the six-member ensem- 
ble as the seasonality of the response (and in particular of 
the decaying phase of EN) is more realistic in these exper- 
iments. However, all of the features described apply to the 
perpetual ENSO experiments as well. 

[ 14 ] Figure 2 shows composites of SSTa, tropopause tem- 
perature, and water vapor, during EPW (top half) and CPW 
(bottom half). In all seasons and for both types of EN, the 
tropopause temperature response resembles a Rossby wave 
over the Indian Ocean and Western Pacific, and a meridional 
dipole straddling the equator over the central and eastern 
Pacific [as in Xie et al, 2012]. However, subtle differences 


exist among the different calendar months and between 
CPW and EPW, and these differences influence the response 
of water vapor to EN. 

[15] We start with a discussion of the seasonality of the 
response during EPW. In NDJF, two factors lead to net cool- 
ing of the CPR at 100 hPa. First, the warming at 100 hPa 
over the warm pool region does not extend as far into the 
Central Pacific as compared to AMJ (Figures 2b and 2d). 
Therefore, only part of the origination region for strato- 
spheric water vapor (i.e., the area enclosed by the green 
isotherm) warms. Second, the CPR is shifted eastward due 
to the warming associated with EPW. EPW induces a net 
warming of the climatological CPR (i.e., the area enclosed 
by a cyan contour), but not of the CPR after it has been 
shifted (i.e., the area enclosed by a green contour). This is 
shown by the dashed red (fixed CPR) and solid red (EN 
CPR) lines in Figure 2i. Finally, the CPR cools in NDJF: the 
temperature of the coldest 10% of the Tropics (i.e., the area 
enclosed by the cyan and green contours) cools by 0.25 K. 

[ 16 ] In contrast, in AMJ, the CPR warms at 100 hPa 
(Figure 2i). First, the warming over the warm pool region 
extends further toward the central Pacific. In addition, dur- 
ing AMJ, the origination region for stratospheric water vapor 
begins to extend toward the Indian subcontinent (cf. green 
isotherm in Figure 2d), and in this region, EPW induces 
warming at 100 hPa. Finally, the CPR warms in AMJ: the 
temperature of the coldest 10% of the Tropics (i.e., the tem- 
perature marked by a green contour) increases by 0.35 K 
during EPW as compared to climatology. In summary, EPW 
causes the CPR to cool in fall and winter, but to warm 
in spring. Hence, we see a transition from a dryer lower 
stratosphere in NDJF to a moister lower stratosphere in 
springtime. 

[17] During CPW, the temperature response at 100 hPa 
is generally similar to that of EPW. However, subtle differ- 
ences exist between the CPW and EPW responses, and these 
differences lead to a dryer springtime stratosphere during 
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Figure 2. Response to EPW and CPW in the ensemble mean of the six-member ensemble forced with observed SSTs only, 
(a-d) EPW and (e-h) CPW. (Figures 2a, 2b, 2e, and 2f) NDJF and (Figures 2c, 2d, 2g, and 2h) AMJ. (left column) Sea 
surface temperature; (middle column) temperature at 100 hPa; (right column) cold point temperature (as in Figure le) and 
zonally averaged water vapor, with contours shown at ±0.1 ppmv, ±0.2 ppmv, ±0.4 ppmv, ±0.6 ppmv, ±0.8 ppmv, and 
±1.5 ppmv. For the middle and right columns, color indicates regions in which the model responses are robust among the 
ensemble members. A dashed red line on Figure 2i shows the average temperature in the region included climatologically 
in the cold point region (the cyan contour; fixed cold point). Figure 2i includes the average cold point temperature if a 
threshold of 6% or 14% is used instead of 10% with a thin line. A solid black contour indicates the climatological position 
of the tropopause in Figures 2j-2k. 


CPW as compared to EPW. First, the magnitude of the tem- 
perature anomalies in the TTL is approximately 50% larger 
during EPW than during CPW in all seasons. The weaken- 
ing of CPW teleconnections relative to EPW teleconnections 
is consistent with the underlying SSTa: the magnitude of 
the SSTa is approximately 50% weaker during CPW than 
during EPW. In addition, the warm anomaly over the warm 
pool region during CPW does not extend as far east into the 
Pacific (this can be seen by comparing the location of the 
node near 180° in Figure 2d versus Figure 2h). The zonal 
extent of the warm anomaly at 100 hPa is related to the 
underlying SST anomalies in each ENSO type, as the posi- 
tive SSTa in the tropical Pacific Ocean during CPW extend 
further to the west than those of EPW (Figures 2c and 2g). 
The upper tropospheric tropical temperature response is also 
westward shifted during CPW as compared to EPW (not 
shown). The net effect is that in AMJ, the CPR warms, 
and the lower stratosphere moistens more during EPW than 
during CPW (cf. Figures 2i and 2k). 


[is] In summary, the net warming of the CPR associated 
with EN peaks in the late spring after the EN event has 
already substantially decayed (Figures le and 2i). Hence, the 
plume of enhanced stratospheric water vapor is delayed until 
the late spring following an EN event. The wanning of the 
CPR is stronger and more pervasive during EPW than during 
CPW (Figures le and 2i). Therefore, the lower stratosphere 
is moister during EPW events. 

4. Conclusions 

[ 19 ] Targeted chemistry-climate model experiments are 
used to demonstrate the importance of seasonality and the 
variant of El Nino event when evaluating the response of 
water vapor to El Nino. There are two main conclusions to 
this study: 

[ 20 ] 1. El Nino events in which sea surface temperature 

anomalies peak in the eastern Pacific lead to significantly 
more water vapor entering the stratosphere as compared 
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Figure 3. The correlation between equatorial (5°S-5°N) 
SST anomalies and zonal-mean water vapor anomalies at 
85 hPA 2 months later. November through February (NDJF) 
SST with January through April (JFMA) water vapor (red), 
and April through June (AMJ) SST with June through 
August (JJA) water vapor (blue). The Nino3 and Nino4 
regions are indicated directly on the figure with cyan and 
green. Correlations whose absolute value exceeds 0.15 are 
statistically significant at the 95% level by a two-tailed 
Student’s t test; these values are indicated on the figure. 


to El Nino events in which the sea surface temperature 
anomaly is found mainly in the central Pacific. The differ- 
ence approaches 0.3 ppmv in certain calendar months. 

[ 21 ] 2. The moistening of the lower stratosphere in 
response to El Nino occurs mainly in spring. In contrast, in 
winter, El Nino leads to a drying of the lower stratosphere. 
The difference between the winter and spring responses 
exceeds 0.5 ppmv. 

[ 22 ] Both of these features are consistent with, and are 
driven by, the temperature anomalies in the Indo-Pacific 
region, and in particular near the cold point, at 100 hPa. 

[ 23 ] The connection between SSTa and stratospheric 
water vapor anomalies is summarized by Figure 3, which 
shows the correlation between stratospheric water vapor and 
tropical SSTa at a lead of 2 months. A 2 month lead is chosen 
as it approximates the timescale for water vapor anomalies 
to reach the 85 hPa level (cf. Figures lc, Id, 2j, and 2k). 
Correlation coefficients exceeding 0.6 are present in the 
tropical Pacific Ocean. During NDJF, La Nina SSTa lead to 
moistening of the stratosphere 2 months later, while during 
AMJ, warming in the eastern Pacific (but not in the cen- 
tral Pacific) leads to moistening of the stratosphere 2 months 
later. For reference, the Nino3 and Nino4 regions are indi- 
cated directly on the figure; these two regions are imperfect 
if one wishes to isolate the SST pattern that most strongly 
leads to moistening of the stratosphere. In addition, Figure 3 
emphasizes the importance of seasonality when considering 
the response of stratospheric water vapor to ENSO. 

[ 24 ] Unfortunately, the available satellite data on strato- 
spheric water vapor is too short to discriminate a statistically 
significant difference between the El Nino variants in obser- 
vations. In addition, the complexity of the sequence of 
physical events leading from SST forcing to stratospheric 
response raises questions about any conclusions based on 
an individual atmospheric GCM. For future work, we plan 
on examining these effects in CCMVal (and Coupled Model 


Intercomparison Project Phase 5) models. Plowever, the 
model output presented here clearly shows that subtle tem- 
perature differences at 100 hPa between winter and spring, 
and between different El Nino variants, can lead to dra- 
matic differences in stratospheric water vapor. As many 
coupled ocean-atmosphere models have trouble reproducing 
the observed diversity of ENSO events in the present-day 
climate [Kim and Yu, 2012], they likely miss these subtle 
differences. If central Pacific El Nino events are becoming 
more frequent (as suggested by Yeh et al. [2009] though dis- 
puted by L ’Heureux et al. [2012]), then the overall impact of 
ENSO on stratospheric water vapor will likely change. 
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